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A  methodology,  based  on  synchrotron  X-ray  nano-tomography  measurements,  is  proposed  to  obtain  3D 
reconstructions  of  porous  supports  for  Solid  Oxide  Fuel  Cell  (SOFC).  The  methodology  has  been  applied  to 
investigate  the  microstructure  of  a  typical  Ni-8YSZ  cermet  substrate  of  an  Anode  Supported  Cell  (ASC). 
Experimental  conditions  have  been  optimised  so  that  a  3D  reconstruction  of  42  |jim  x  42  |jim  x  105  pirn 
has  been  computed  with  a  voxel  size  of  60  nm. 

The  3D  reconstruction  has  been  numerically  analysed  showing  a  quasi  isotropy  over  the  medium.  A 
special  attention  has  been  paid  to  compute  the  morphological  properties  controlling  the  gas  diffusion 
through  the  support.  In  this  frame,  it  has  been  found  that  a  volume  as  large  as  35  p,m  x  35  (Jim  x  35  (Jim  is 
required  to  be  statistically  representative  for  the  whole  substrate.  More  particularly,  the  tortuosity  factor 
r  has  been  calculated  on  the  basis  of  finite  element  simulations.  The  effect  of  the  boundary  conditions 
taken  for  these  simulations  has  been  investigated.  In  the  direction  perpendicular  to  the  anode/electrolyte 
interface,  the  tortuosity  factor  has  been  determined  to  rz  =  2.8+q \. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  Oxide  Fuel  Cells  (SOFCs)  are  electrochemical  devices  that 
ensure  the  conversion  of  fuel  chemical  energy  into  electricity  at 
high  temperature.  Since  they  present  high  efficiencies  and  large  fuel 
flexibility,  SOFCs  appear  as  an  attractive  technology  for  economic 
power  generation.  The  active  part  of  SOFC  corresponds  to  a  three¬ 
layered  structure,  constituted  by  a  dense  electrolyte  sandwiched 
between  two  porous  electrodes.  The  most  common  materials  used 
for  this  membrane  are  Lanthanum  Strontium  Manganite  (LSM)  for 
the  cathode,  Yttria  Stabilized  Zirconia  (YSZ)  for  the  electrolyte,  and 
Nickel/YSZ  cermet  (Ni-YSZ)  for  the  anode. 

Nowadays,  efforts  are  ongoing  to  operate  the  SOFCs  at  interme¬ 
diate  temperature  (T <  800  °C)  in  order  to  reduce  the  material  costs 
and  simultaneously  improve  the  cell  reliability.  Anode  Supported 
Cells  (ASCs)  have  been  designed  for  such  application  [1,2].  The  ASC 
configuration  is  based  on  a  thick  highly  porous  anodic  substrate, 
on  which  a  thin  electrolyte  and  the  functional  electrodes  are  lay¬ 
ered.  The  substrate,  which  is  used  as  the  cell  structural  support, 
presents  a  coarse  microstructure  to  facilitate  the  gas  transporta¬ 
tion;  whereas  the  fine  microstructure  of  the  functional  layer  is 
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optimised  for  the  electrochemical  reactions  [3,4].  Thanks  to  its  con¬ 
figuration,  ASC  allows  reducing  the  ohmic  losses  in  the  electrolyte 
and  exhibits  high  power  densities.  Furthermore,  it  has  been  shown 
that  this  design  is  well  adapted  to  operate  with  natural  gas  [5]. 

Flowever,  the  microstructure  features  of  the  electrodes  are  key 
parameters  which  determine  the  cell  polarization  resistance,  and 
hence,  the  electrochemical  performances  [6-8].  The  evolution  of 
these  microstructure  parameters  after  a  long  term  operation  or 
after  redox  cycling  could  also  explain  the  degradation  in  electro¬ 
chemical  performances  [9-12]. 

More  particularly,  it  has  been  shown  that  gas  diffusion  into  the 
ASC  substrate  can  induce  high  concentration  overpotentials,  which 
can  limit  the  cell  efficiency  for  operation  at  elevated  current  density 
[13,14].  It  is  worth  noting  that  the  gas  transport  up  to  the  electro¬ 
chemical  active  sites  is  controlled  by  the  morphological  properties 
of  the  electrode  supports  [7,15],  more  particularly  by  the  3D  topol¬ 
ogy  of  the  pores  network.  These  parameters  also  affect  the  methane 
conversion  rate  for  SOFC  operating  under  direct  internal  reforming 
[16].  As  a  consequence,  the  knowledge  of  these  structural  prop¬ 
erties  is  essential  to  predict  the  mass  transport  into  the  porous 
substrate  which  is  linked  to  the  ASC  efficiency. 

For  all  these  reasons,  several  efforts  have  been  recently  paid  (i) 
to  obtain  3D  reconstructions  of  real  electrodes  and  (ii)  to  charac¬ 
terise  their  morphological  and  fluid  properties  [  1 1,17-29].  Focused 
Ion  Beam  coupled  with  a  Scanning  Electron  Microscope  (FIB-SEM) 
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[11,17-25]  and  X-ray  nano-tomography  [26-29]  methods  have 
been  widely  employed  for  this  purpose.  It  has  been  demonstrated 
that  both  techniques,  radically  different  in  their  basic  principles, 
stand  to  obtain  3D  reconstruction  of  relatively  small  volume  (typ¬ 
ically  close  to  10  pan  x  10  pun  x  10  pan  [23])  at  high  resolution 
(as  low  as  lOnm  for  FIB-SEM  technique  [11]  and  45  nm  for  X- 
ray  tomography  [29]).  They  have  been  successfully  applied  to 
determine  the  microstructure  properties  of  functional  electrodes. 
Because  of  the  fine  microstructure  of  such  layers,  it  has  been  shown 
that  the  volumes  of  the  reconstructions  are  sufficiently  large  to 
reach  a  Representative  Volume  Element  (RVE)  of  the  analysed 
medium  [18-20].  For  instance,  Wilson  et  al.  [18]  have  stated  that 
the  analysis  of  volumes  higher  than  >500  |xm3  are  sufficient  to 
provide  microstructure  properties  that  are  representative  for  the 
overall  electrode.  This  result  is  confirmed  by  Cronin  et  al.  [19]  who 
suggested  that  volumes  of  ~1000fxm3  are  large  enough  to  pro¬ 
vide  good  statistics.  These  volumes  appear  to  be  slightly  lower 
than  the  one  recommended  by  Metcalfe  et  al.  [30]:  by  a  numer¬ 
ical  approach,  these  authors  have  established  that  the  size  of  the 
domain  should  be  14  times  the  mean  particles  diameter  to  ade¬ 
quately  represent  the  electrode  structure.  It  can  be  noticed  that  Joos 
et  al.  [24]  have  succeeded  to  reconstruct  the  entire  cathode  thick¬ 
ness  of  an  ASC  (volume  of  30.7  [xm  x  5.25  |xm  x  1 5.75  |xm  obtained 
by  FIB-SEM).  The  largest  reconstruction  was  achieved  by  Kanno 
et  al.  [25]  for  a  conventional  Ni-YSZ  anode  deposited  onto  a  thick 
electrolyte:  the  authors  have  investigated  a  reconstructed  struc¬ 
ture  as  large  as  45.8  |xm  x  14.5  |xm  x  26.2  |xm.  However,  at  this 
time,  no  study  has  been  devoted  to  the  ASC  substrate,  probably 
because  its  coarse  microstructure  requires  an  analysis  over  a  larger 
volume. 

The  present  study  aims  at  overcoming  this  limitation  and  is  ded¬ 
icated  to  the  reconstruction  of  a  standard  Ni-ZSZ  cermet  support. 
For  this  purpose,  hard  X-ray  nano-tomography  has  been  used  to 
obtain  valuable  3D  imaging  of  the  microstructure  at  a  relevant  scale. 
The  3D  reconstruction  has  been  analysed  by  using  computational 
modelling  in  order  to  determine  the  morphological  properties 
controlling  the  gas  diffusion  in  the  substrate.  These  parameters  cor¬ 
respond  to  the  fluid  contiguity,  the  porosity,  the  mean  pore  radius 
and  the  tortuosity  factor. 


2.  Experimental 

2. 1.  Materials  and  sample  preparation 

Commercial  ASC  supplied  by  the  Forschungszentrum  Jiilich  (FZJ) 
research  centre  [4]  has  been  used  for  the  experiments.  The  anode 
includes  a  thin  Anode  Functional  Layer  (AFL)  deposited  on  a  thick 
Anode  Current  Collector  (ACC).  The  ACC  is  ~1  mm  thick  and  is  used 
as  the  structural  support  of  the  cell.  As  received,  the  anode  is  consti¬ 
tuted  of  56  wt.%  NiO  and  44  wt.%  8YSZ  (Zr02  stabilized  with  8  mol% 
Y203 ).  The  AFL  and  the  electrolyte  (8YSZ)  are  thin  layers  of  ~1 0  |xm 
whereas  the  cathode  thickness  is  ~50  [xm. 

A  narrow  strip  of  cell  has  been  prepared  by  laser  cutting.  This 
strip  was  polished  up  to  achieve  a  thin  plate  of  100  fxm  in  thick¬ 
ness.  This  plate  was  finally  sliced  along  the  cell  thickness  by  wire 
cutting  to  obtain  the  multi-layered  samples  for  the  X-ray  tomog¬ 
raphy  experiments.  The  geometry  and  dimensions  of  the  samples 
are  depicted  in  Fig.  1 :  they  correspond  to  “needles”  whose  the  rev¬ 
olution  axis  is  perpendicular  to  the  cell  layers. 

Before  the  tomography  experiments,  the  samples  were  heated 
up  to  T=  800  °C  at  a  rate  of  5  °C  min-1  under  (2%H2  +  Ar).  They  were 
maintained  in  these  conditions  during  48  h  in  order  to  achieve  the 
complete  reduction  of  NiO  into  Ni.  This  procedure  has  been  val¬ 
idated  by  X-ray  diffraction  indicating  that  all  the  Ni  particles  are 
present  under  metallic  form  after  the  reduction  treatment. 


Z  (revolution  axis) 


Fig.  1.  Schematic  representation  of  the  samples  prepared  for  the  X-ray  tomography 
experiments.  The  zoom  at  one  extremity  of  the  sample  corresponds  to  a  Scanning 
Electron  Microscope  (SEM)  image. 

The  porosity  of  the  cermet  support  was  measured  by 
Archimedes’  method.  The  porosity,  which  is  denoted  e  and  defined 
as  the  fraction  of  pores  volume,  was  determined  to  s  ~  47%  for  the 
material  in  its  Ni-8YSZ  reduced  state.  The  specific  surface  area  was 
also  measured  by  the  BET  method  [31  ].  For  this  purpose,  pieces  of 
cermet  support  («7  mm  x  5  mm  x  1  mm)  were  prepared  and  char¬ 
acterised  in  an  Analyser  SA  3100  Beckman  Coulter-NG  apparatus. 
The  specific  surface  area,  denoted  Sp  and  defined  as  the  surface  area 
normalised  by  the  pore  volume,  was  determined  for  the  Ni-8YSZ 
support  to  Sp  =  2.51  |xm-1 . 

2.2.  X-ray  tomography  procedure 

To  produce  large  3D  reconstruction  of  SOFCs,  hard  X-rays  must 
be  used  to  prevent  from  excessive  absorption  of  the  beam  within 
the  sample.  Given  the  critical  dimensions  of  SOFCs,  high  resolu¬ 
tion  data  are  required,  fostering  the  use  of  brilliant  hard  X-ray 
sources  as  provided  by  3rd  generation  synchrotrons.  Image  acqui¬ 
sition  has  been  carried  out  at  the  European  Synchrotron  Radiation 
Facility  (ESRF)  beamline  ID22NI,  based  on  a  holotomographic  scan¬ 
ning  scheme  [32].  The  beam  energy  is  fixed  to  29keV  with  1.5% 
bandwidth.  Such  high  energy  is  required  to  obtain  a  sufficient 
transmitted  X-ray  flux  to  investigate  a  large  volume  of  material. 
Furthermore,  this  energy  is  high  enough  to  get  ideal  phase  retrieval 
conditions.  Indeed,  scans  at  1 7  keV  confirmed  that,  with  our  “thick” 
sample  of  Ni-8YSZ,  the  final  image  quality  was  not  good  enough 
due  to  the  strong  absorption  of  the  X-ray  beam  and  also  because  of 
the  strong  contrast  in  density  of  such  porous  material.  The  29  keV 
X-ray  beam  is  focused  down  to  a  spot  below  lOOnm  in  horizontal 
and  vertical  directions  using  Kirkpatrick-Baez  mirrors.  The  beam 
then  diverges  to  illuminate  the  sample  placed  downstream  from 
the  focal  spot.  A  2D  X-ray  detector  is  placed  further  downstream  to 
record  a  magnified  image  of  the  sample  [33,34].  In  this  procedure, 
images  are  not  standard  radiographs  but  rather  holograms:  4  dif¬ 
ferent  tomographic  scans  at  4  sample-to-detector  distances,  each 
scan  consisting  of  a  series  of  2000  projections  over  a  turn,  must  be 
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Fig.  2.  Representation  of  the  three-dimensional  reconstructed  microstructure  of  the 
studied  Ni-8YSZ  support. 
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Fig.  3.  Porosity  s  plotted  as  a  function  of  the  cube  size  i  extracted  from  the  whole 
reconstruction. 


This  section  of  the  paper  is  focused  on  the  determination  of  these 
three  structural  parameters  from  the  3D  reconstruction. 

3.1.  Fluid  contiguity  and  porosity 


acquired.  The  series  of  4  tomographic  scans  serves  as  an  input  to  a 
phase  retrieval  procedure  that  produces  projection  images  whose 
gray  levels  are  proportional  to  the  projected  electron  density.  The 
projected  electron  density  images  are  then  used  as  an  input  to  a 
3D  volume  reconstruction  algorithm.  The  full  sample  imaging  lasts 
about  90  min.  The  voxel  size  used  is  60  nm,  giving  the  field  of  view 
of  about  120  pirn,  allowing  analysing  samples  with  cross-sections 
close  to  this  value. 

2.3.  Three-dimensional  reconstruction 

In  order  to  address  the  full  microstructure  of  the  Ni-8YSZ  cer¬ 
met,  it  would  be  required  to  separate  the  gas  phase  as  well  as 
the  metallic  and  the  ceramic  phases  of  the  composite.  By  X-ray 
nanotomography,  Grew  et  al.  [26]  and  Shearing  et  al.  [28]  have 
shown  that  this  can  be  achieved  by  using  differential  absorption 
imaging  across  the  Ni  K-absorption  edge  (8.333  keV).  However, 
this  low  beam  energy  is  detrimental  to  analyse  large  volume 
due  the  strong  absorption  of  the  materials.  For  instance,  Nelson 
et  al.  [35]  have  used  cylindrical  sample  approximately  lOpim  in 
diameter  and  10-20  pun  in  height.  As  mentioned  previously,  such 
volumes  are  expected  to  be  not  large  enough  regarding  to  the 
coarse  microstructure  of  SOFC  support.  Higher  beam  energy  is  then 
required.  However,  under  this  condition,  Ni  and  8YSZ  exhibits  sim¬ 
ilar  X-ray  absorption  properties  which  prevent  from  solid  phase 
separation. 

The  3D  reconstruction  of  the  microstructure  obtained  in  this 
work  is  illustrated  in  Fig.  2.  A  volume  of  42  pum  x  42  pum  x  105  pum 
is  obtained  in  which  the  solid  and  pore  phases  are  identified.  This 
reconstructed  volume  is  located  within  the  ACC  just  above  the  func¬ 
tional  layer. 

3.  Numerical  analysis:  quantification  of  microstructure 
parameters 


The  connectivity  of  the  gas  phase  can  be  represented  by  the  fluid 
contiguity  which  is  defined  by  the  volumetric  percentage  of  the 
connected  pores  [26].  An  in-house  program  has  been  implemented 
in  the  commercial  Matlab®  software  to  calculate  this  parameter. 
The  program  has  been  developed  by  fulfilling  the  following  rules: 
two  voxels  of  pores,  that  shared  one  common  face,  were  considered 
as  connected  to  form  a  cluster;  whereas  two  voxels  only  linked  by 
an  edge  or  a  vertex  were  considered  as  disconnected.  As  described 
in  Ref.  [26],  the  percolating  volume  of  pores  was  then  identified  by 
removing  from  the  network  the  isolated  clusters  which  are  either 
embedded  as  closed  porosity  in  the  bulk  or  emerge  only  at  one  edge 
of  the  3D  reconstruction.  The  3D  pores  network  of  the  reduced  cer¬ 
met  has  been  characterised  by  this  method.  The  numerical  analysis 
has  revealed  that  96%  of  the  pore  phase  is  percolating.  This  result 
means  that  the  main  part  of  the  pores  is  connected  in  one  single 
network,  suggesting  that  the  porous  medium  exhibits  a  good  ability 
regarding  the  gas  flow. 

Porosity  has  been  also  determined  for  a  set  of  sub-volumes, 
which  have  been  randomly  extracted  from  the  whole  3D  recon¬ 
struction.  These  sub-volumes  correspond  to  cubes  for  which  the 
edge  length  i  ranges  from  0.3  pum  to  40.8  pum.  Considering  domains 
with  the  same  size  i ,  a  sufficient  large  number  of  analyses  allows 
determining  the  average  porosity  and  its  standard  deviation.  The 
mean  porosity  is  given  in  Fig.  3  and  agrees  the  value  provided  by 
the  Archimedes’  method  (i.e.  47%:  see  Section  2.1 ).  This  result  tends 
to  confirm  the  accuracy  of  the  reconstruction. 

As  illustrated  in  Fig.  3,  the  standard  deviation  from  the 
mean  porosity  value  decreases  with  increasing  the  volume 
size.  The  dispersion  becomes  negligible  for  volume  higher  than 
35  pum  x  35  pum  x  35  pum.  In  other  words,  this  volume  is  suffi¬ 
ciently  large  to  be  statistically  representative  of  the  heterogeneous 
material  and  can  be  taken  as  the  threshold  corresponding  to  the 
Representative  Volume  Element  (RVE)  for  the  investigated  prop¬ 
erty. 


The  multi-component  transport  in  the  porous  cermet  structure 
is  usually  modelled  in  the  frame  of  the  Mean  Transport  Pore  Model 
(MTPM)  [15,16].  This  model  is  based  on  the  assumption  that  the 
porous  medium  can  be  described  by  three  microstructure  param¬ 
eters:  the  porosity,  the  tortuosity  factor  and  the  mean  pore  radius. 


3.2.  Surface  area  and  mean  pore  radius 

The  specific  surface  area  is  defined  as  the  surface  area  nor¬ 
malised  by  the  phase  volume  [18].  This  parameter  has  been 
calculated  for  the  porosity.  For  this  purpose,  each  voxel  face  of  the 
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Fig.  4.  Covariograms  calculated  from  the  3D  reconstruction.  The  covariance  K  has 
been  computed  in  the  x,  y,  z  directions. 
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3D  reconstruction  has  been  examined:  interfaces  between  solid  and 
gas  have  been  identified  and  then  enumerated  in  order  to  calculate 
the  fluid/solid  area  thanks  to  the  surface  of  the  voxel  face  (i.e.  i 2oxd ). 
It  is  worth  noting  that  the  value  calculated  by  this  method  requires 
a  rescaling  in  order  to  estimate  the  real  surface.  Indeed,  as  the  struc¬ 
ture  is  described  by  cubic  elements  (i.e.  the  voxels),  the  surface  of 
objects  included  in  the  digitized  reconstruction  is  overestimated 
[36].  Under  simple  geometrical  considerations,  it  can  be  shown  that 
the  scaling  corresponds  to  a  factor  of  tt/6,  if  the  objects  are  spheri¬ 
cal  or  randomly  oriented  (indeed,  in  this  case,  the  factor  is  equal  to 
the  ratio  between  the  surfaces  of  a  sphere  and  a  cube).  As  a  conse¬ 
quence,  it  can  be  reasonably  supposed  the  real  specific  pore  surface 
area  ranges  between  the  rough  data  deduced  directly  from  the 
numerical  reconstruction  and  the  processed  value  obtained  after 
the  scaling  by  a  factor  tt/6:  1.80  [xm-1  <SP<3.45  [xnrr1.  It  worth 
noting  that  the  median  value,  i.e.  2.62  |xm-1,  is  found  to  be  in 
good  agreement  with  the  data  estimated  by  the  BET  method,  i.e. 
2.51  |xm-1  (see  Section  2.1). 

It  is  assumed  that  an  assemblage  of  spherical  pores  (all  with  the 
same  size)  is  representative  for  the  material  porosity.  In  this  case, 
the  specific  surface  area  can  be  linked  to  the  mean  pore  radius  Rp 
according  to  [37]: 

s>=h  (1) 

This  relation  leads  to  determine  an  average  pore  diameter  of 
2  x  Rp  =  2.4  |x  m.  In  order  to  check  this  value,  the  covariogram  I<(h) 
of  the  reconstruction  has  been  computed  from  the  reconstruction 
[38].  It  can  be  expressed  as  a  measure  of  the  intersection  between 
the  porosity  of  the  material  and  its  image  after  a  translation  h  of  the 
medium.  Considering  the  x-direction,  I< can  be  written  as  followed: 


I<(h) 

with 


k{x,  y,  z )  x  k((x  +  h),  y,  z )  x  dV , 

J  1  if  x  belongs  to  the  pore  phase 
1  0  if  x  is  in  the  solid  phase 


(2) 


Based  on  the  3D  tomographic  images,  covariograms  have  been 
computed  in  the  three  directions  and  are  provided  in  Fig.  4.  As 
expected,  I<(h)  is  equal  to  the  value  of  the  porosity  for  h  =  0  and  tends 
to  an  asymptotic  theoretical  value  corresponding  to  the  square  of 
the  porosity:  lim I<(h)  =  e2  [38].  The  distance  h*  to  reach  this  limit 

h — ^  oo 

corresponds  to  the  range  of  the  covariance  and  can  be  ascribed  to 
the  characteristic  length  of  the  porosity  [38].  A  criterion  chosen  as 
((K(h*)  -  s2)ls2)<  0.05  leads  to  an  estimation  of  the  characteristic 


Fig.  5.  Schematic  representation  of  sinuous  but  parallel  pores  to  define  the  tortu¬ 
osity. 


size  for  the  pores  of  h*  ~  3  |xm.  This  value  is  consistent  with  the 
mean  pore  diameter  2 Rp  assessed  previously  by  assuming  a  pack¬ 
ing  of  spherical  pores.  As  a  consequence,  it  can  be  concluded  that 
the  mean  pore  radius  of  the  analysed  substrate  ranges  between  1.2 
and  1.5  |xm. 

3.3.  Tortuosity  factor 

Definition  and  numerical  method  -  Tortuosity  f  can  be  calculated 
for  the  three  directions  of  the  3D  reconstruction.  For  the  sake  of 
simplicity,  this  parameter  can  be  easily  expressed  considering  a  2D 
representation  as  depicted  in  Fig.  5.  In  this  case,  the  tortuosity  is 
defined  as  the  ratio  between  the  effective  path  length  ^ eff  and  the 
length  of  the  sample  lx  (Fig.  5): 


Considering  a  porous  medium  as  a  bundle  of  sinuous  but  parallel 
capillaries,  it  can  be  demonstrated  that  the  fluid  effective  diffusion 
coefficient  Df ^  is  given  by  [39]: 

of  =  ^fbutk  with  Xi  = 

i  =  bundle-direction  (x,y  or  z)  (4) 


The  term  r2-  is  usually  named  as  the  tortuosity  factor  and  Dbuik 
denotes  the  diffusion  coefficient  of  the  gas.  In  porous  materials  such 
as  SOFC  electrodes,  there  are  many  pathways  that  are  entangled 
together.  For  this  reason,  it  is  impossible  and  not  relevant  to  define 
an  effective  transport  path  length,  ^eff*  and  hence,  a  “geometrical” 
tortuosity  f  of  the  medium  [22,39].  In  other  terms,  the  effective 
diffusion  coefficient  cannot  be  determined  from  geometrical  con¬ 
siderations.  However,  provided  that  Df ^  can  be  calculated  by  other 
methods,  it  becomes  possible  to  calculate  the  empirical  ratio  of  the 
porosity  to  a  “reliable  tortuosity  factor”,  (e/r,*)  =  (DeP / Dbuik ),  which 
takes  into  account  the  effects  of  structure. 

Several  methods  have  been  proposed  to  evaluate  the  tortuosity 
factor  through  the  calculation  of  an  effective  diffusion  coefficient. 
The  tortuosity  factor  can  be  calculated  by  using  the  Lattice  Boltz- 
man  Method  (LBM)  [20,26].  Iwai  et  al.  [20,22]  have  statistically 
calculated  this  parameter  from  the  random  walk  process.  Wilson 
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Fig.  6.  Concentration  field  in  the  percolated  porosity  (volume  of 
36  (Jim  x  36  (Jim  x  36  jxm  extracted  from  the  3D  reconstruction,  two  concen¬ 
trations  imposed  at  z  =  0  and  z  =  36  |xm  as  boundary  condition,  Element  (or  voxel) 
size  lvoxei  =  240  nm,  2.2  x  106  nodes). 


et  al.  [17]  and  Joos  et  al.  [24,40]  have  solved  the  Laplace’s  equation 
for  diffusion  within  the  microstructure  by  using  the  Finite  Element 
Method  (FEM).  Although  the  boundary  conditions  chosen  for  these 
simulations  can  introduce  a  bias  in  the  assessment  of  the  effective 
property  [38],  this  effect  has  not  been  estimated  for  SOFC  electrode 
yet. 

In  this  work,  mass  transfer  by  diffusion  in  the  gas  phase  of  the 
microstructure  has  been  simulated  by  FEM.  The  Cast3M  solver  [41  ] 
has  been  used  for  this  purpose.  The  mesh  has  been  built  by  convert¬ 
ing  each  voxel  of  the  open  porosity  into  an  8-nodes  cubic  element. 
Thanks  to  the  methodology  used  to  identify  the  percolated  gas 
phase  (see  Section  3.1),  no  element  can  be  connected  to  a  neigh¬ 
bour  by  an  edge  or  a  vertex.  This  method  prevents  from  reducing 
locally  to  zero  the  size  of  the  section  for  gas  transport  which  would 
act  as  a  geometrical  singularity  in  the  mesh. 

Only  cubic  volumes  extracted  from  the  reconstruction  have 
been  analysed  in  order  to  avoid  any  anisotropic  dimensional  effect 
on  the  result.  Simulations  have  been  carried  out  by  imposing  a  uni¬ 
form  concentration  or  molar  flux  onto  two  opposite  faces  of  the 
domain  ( i.e .  the  opposite  external  faces  of  the  cube  perpendicular 
to  the  three  x,y  or  z  axes).  In  the  other  part  of  the  mesh,  no  boundary 
condition  has  been  applied.  It  can  be  noticed  that  this  condition  is 
equivalent  to  obtain  a  zero  molar  flux  through  the  solid/gas  inter¬ 
face  in  the  real  sample.  Over  the  modelled  structure,  the  Laplace’s 
equation  for  diffusion  is  solved  by  the  FEM: 

f  d2C  d2C  d2C  1 

div[-Dbuik  X  grade]  =  0«  -Dbulk  j_  +  _  +  _j=0  (5) 

where  C  is  the  concentration  at  the  location  (x,  y,  z)  and  Dbuik  is 
supposed  independent  of  C.  As  an  illustration,  a  simulated  con¬ 
centration  field  in  the  percolated  porosity  is  plotted  in  Fig.  6 
for  a  volume  of  36  [xm  x  36  |xm  x  36  |xm  extracted  of  the  cermet 
reconstruction.  It  can  be  noticed  that  the  mesh  for  this  simulation 
contains  2.2  x  106  nodes. 

Once  the  concentration  field  is  known,  the  effective  diffusion 
coefficient  can  be  easily  deduced  from  the  net  flux  N;  entering  or 


exiting  from  the  domain: 

Nt  _  Deff  (Q=a/>  -  (Cm,) 

AjAk  '  A  i 


1 

AjAk 


l 


Q  x  ds 


(i  =  x,y 


with  (Q) 

or  z;  j  ±  i  and  k  j=  i,j) 


(6) 


where  Aj=Ak  =  Ai  =  £  is  the  edge  length  of  the  simulated  cubic 
volume  and  S  its  surface  (S  =  £  2 ).  The  tortuosity  factor  is  then  deter¬ 
mined  by  combining  Eqs.  (4)  and  (6). 

Sensitivity  analysis  -  A  sensitivity  analysis  has  been  performed 
on  the  tortuosity  factor  by  changing  either  the  voxel  size  lvoxei 
or  the  simulated  volume  length  i.  A  specific  procedure  has  been 
implemented  in  Matlab®  to  decrease  the  resolution  of  the  3D  recon¬ 
struction.  In  the  method,  a  cube  of  n3  (n  =  2,  3,  4, . . .)  initial  voxels 
have  been  merged  together  to  form  a  new  cubic  element.  In  the  case 
where  such  a  coarse  voxel  contained  sub-voxels  belonging  to  the 
solid  and  pore  phases,  the  nature  of  the  new  element  was  deter¬ 
mined  as  followed:  it  was  attributed  to  the  gas  phase  when  the 
pore  sub-voxels  were  predominant  within  the  initial  set  of  sub¬ 
voxels.  Conversely,  when  they  were  found  in  minority,  the  new 
element  was  considered  belonging  to  the  solid  phase.  In  case  of 
strict  equality,  the  nature  of  the  element  was  chosen  randomly. 

The  voxel  size  dependence  has  been  investigated  on  a  volume 
of  9  |xm  x  9  |xm  x  9  |xm  extracted  from  the  3D  reconstruction.  The 
simulations  have  been  carried  out  to  calculate  the  tortuosity  fac¬ 
tor  rx  according  to  the  x-direction.  Two  uniform  concentrations 
have  been  imposed  on  the  left  and  right  surfaces  of  the  mod¬ 
elled  volume  respectively  at  x  =  0|xm  and  x  =  9|xm.  In  Fig.  7,  the 
error  on  the  tortuosity  factor  has  been  plotted  as  a  function  of 
the  voxel  length.  It  can  be  seen  that  the  size  of  the  discretization 
must  be  lower  than  lvoxei  <  300  nm  to  avoid  any  error  higher  than 
100  x  |rx  -  r^°xef=60nm|  /r ^oxe/=60nm  >  4.1%.  It  has  been  checked 
that  this  result  remains  valid  for  the  tortuosity  factors,  xy  and  rz, 
determined  according  to  the  y-  and  z-direction.  It  is  worth  noting 
that  this  criterion  will  be  fulfilled  for  all  the  simulations  presented 
in  subsequent  part  of  this  article. 

The  tortuosity  factors  have  been  also  calculated  by  increas¬ 
ing  the  simulated  domain  from  9  |xm  x  9  |xm  x  9  |xm  up  to 
41 .4  [xm  x  41 .4  |xm  x  41 .4  |xm.  In  this  studied  case,  it  can  be  noticed 
the  porosity  decreases  with  increasing  investigated  volumes  from 
~55%  down  to  ~47%.  Simulations  have  been  performed  by  impos¬ 
ing  two  uniform  concentrations  at  the  boundaries.  Results  are 
plotted  in  Fig.  8.  It  is  found  that  the  values  become  steady  for 
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Domain  size  f  (jam) 

Fig.  8.  Tortuosity  factors  xx,  xy,  xz  plotted  as  a  function  of  the  modelled  volume 
length  l.  The  label  gives  the  size  of  the  voxels,  lvoxei,  for  which  the  mesh  has  been 
built.  Simulations  have  been  performed  for  two  imposed  concentrations  at  boundary 
(i.e.  the  opposite  surface  of  the  domain). 


Table  1 

Tortuosity  factors  calculated  for  three  disconnected  volumes  of 
36  |jim  x  36  p,m  x  36  p,m  extracted  from  the  3D  reconstruction  (uniform 
concentration  imposed  at  boundary,  lvoxei  =  240  nm). 


Zone  no.  1 

Zone  no.  2 

Zone  no.  3 

Mean  value 

rx 

1.9 

2.5 

2.5 

2.3 

Xy 

2.9 

2.8 

2.6 

2.8 

Tz 

2.6 

2.7 

2.6 

2.6 

volumes  larger  than  35  |jim  x  35  |jim  x  35  |jim.  This  result  means 
that  the  size  of  the  RVE  relevant  for  the  tortuosity  factor  is  close 
to  the  one  required  for  the  porosity.  As  a  consequence,  this  volume 
can  be  considered  as  a  threshold  from  which  the  reconstruction 
becomes  statistically  representative  for  the  porous  cermet  support. 

Spatial  averaging  -  The  tortuosity  factors  rx ,  ry  and  rz  have  been 
computed  for  three  volumes  of  36  pan  x  36  pan  x  36  pm  taken  at 
different  location  along  the  height  of  the  3D  reconstruction  (i.e. 
according  to  the  vertical  z-axis).  It  is  worth  noting  that  the  extracted 
domains  are  disconnected  from  one  to  each  other.  The  simulations 
have  carried  out  by  imposing  two  uniform  concentrations  at  the 
boundaries. 

The  values  of  the  tortuosity  factors  in  each  direction  are  reported 
in  Table  1 .  It  can  be  noticed  that  the  scattering  on  rx,  ry  or  rz  values 
is  narrow,  i.e.  the  tortuosity  factors  do  not  depend  on  the  location 
of  the  analysed  volume  in  the  support.  This  result  confirms  that 


volumes  as  large  as  35  pan  x  35  pan  x  35  pm  are  required  to  deter¬ 
mine  the  effective  transport  properties  of  the  Ni-YSZ  substrate. 

Moreover,  a  slight  difference  is  noticeable  between  rx,  ry  and  rz 
with  a  similar  trend  for  the  three  modelled  volumes:  ry  >  rz  >  rx.  It 
can  be  inferred  from  this  result  that  the  medium  presents  a  slight 
anisotropy  regarding  to  the  diffusion  of  gaz. 

Effect  of  boundary  conditions  -  The  impact  of  boundary  condi¬ 
tions  has  been  investigated  for  the  three  following  cases  depicted 
in  Table  2: 

(i)  The  first  case  has  been  used  in  the  previous  section:  two  uni¬ 
form  concentrations  are  imposed  on  the  opposite  surfaces  of 
the  volume  (at  z  =  0  and  z  =  i ). 

(ii)  The  second  case  corresponds  to  an  intermediate  boundary  con¬ 
dition  where  one  uniform  concentration  is  imposed  on  the 
lower  surface  (at  z  =  0)  and  a  flux  is  applied  on  the  upper  surface 
(at  z=i). 

(iii)  In  the  last  condition,  two  fluxes  are  applied  on  the  lower  and 
upper  surfaces  of  the  analysed  cube.  It  can  be  noticed  that,  in 
this  case,  the  concentration  on  one  point  of  the  structure  needs 
to  be  defined. 

The  results  on  rz  are  summarized  in  Table  2.  Although  the 
simulated  volume  is  slightly  higher  than  the  estimated  RVE  (i.e. 
35  pm  x  35  pm  x  35  pm),  the  three  investigated  cases  lead  to  a 
significant  difference  on  the  tortuosity  factor.  This  phenomenon 
is  caused  by  some  edged  effects  which  arise  when  the  boundary 
conditions  are  applied  on  a  limited  domain  [38].  From  a  theoreti¬ 
cal  point  of  view,  these  effects  must  vanish  for  an  infinite  volume, 
where  surfaces  of  isoconcentration  values  can  be  considered  to  be 
plane  and  perpendicular  to  the  cube  edges.  A  close  examination 
of  Fig.  6  shows  that  it  is  not  the  case  at  the  scale  of  the  analysed 
volume.  In  practice,  the  first  and  the  last  cases  give  respectively 
a  lower  and  upper  bound  for  the  tortuosity  factor.  Furthermore, 
it  can  be  noticed  that  the  intermediate  case  corresponds  to  the 
closest  conditions  regarding  to  the  SOFC  application.  Indeed,  the 
hydrogen  electro-oxidation  controls  the  H2  and  H20  fluxes  at  the 
anode/electrolyte  interface;  whereas  the  fuel  flowing  along  the 
gas  distributor  imposes  the  concentrations  at  the  other  side  of  the 
cermet.  As  a  consequence,  it  can  be  proposed  to  use  in  SOFC  mod¬ 
els  a  tortuosity  factor  of  rz  =  2.8  for  the  cermet  support  (with  a 
upper  bound  of  t™ax  =  3.1  and  a  lower  bound  of  r™in  =  2.6  denoted 
rz  =  2.8+“). 

4.  Discussion 

RVE  size  and  voxel  resolution  -  Since  they  control  the  diffusion  in 
the  thick  substrate,  the  morphological  properties  of  SOFC  supports 
constitute  key  parameters  regarding  to  the  optimisation  of  SOFCs 
performance.  A  methodology  is  proposed  in  this  paper  to  assess 


Table  2 

Tortuosity  factors  calculated  considering  three  different  boundary  conditions  (volumes  of  36  jxm  x  36  |xm  x  36  jxm  and  Cvoxei  =  240  nm). 


(i)  Two  concentrations  imposed  atz  =  0  and  (ii)  One  concentration  imposed  atz  =  0  and  one  flux  (iii)  Two  fluxes  imposed  atz  =  0  and 

z  =  36  pan  imposed  atz  =  36  p,m  z  =  36  pan 


xz  2.6 


2.8 


3.1 
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these  parameters  in  the  frame  of  a  3D  reconstruction  obtained 
from  nano-tomography  measurements.  Applied  on  a  commercial 
Ni-YSZ  support,  it  has  been  shown  that  bigger  volumes  than  the 
ones  considering  for  the  functional  electrodes  need  to  be  ana¬ 
lysed.  Indeed,  reconstruction  as  large  as  35  pan  x  35  pan  x  35  pan 
is  required  to  be  statistically  representative  of  the  whole  substrate 
regarding  to  the  fluid  properties.  Moreover,  it  is  found  that  the 
ratio  between  the  size  of  the  RVE  and  the  mean  pore  diameter  is 
equal  to  (£^(=35  pm)/(2Rp  )  (=  2.4  p,m))  =  14.6.  This  result  is  con¬ 
sistent  with  the  conclusion  established  by  Metcalfe  et  al.  [30]  who 
have  found  on  the  basis  of  a  purely  numerical  study  that  “a  cubic 
envelope  having  a  minimum  length  of  14  times  the  mean  parti¬ 
cle  diameter  was  necessary  to  adequately  represent  the  electrode 
structure”. 

Conversely  to  the  size  of  the  RVE,  the  resolution  of  the  3D  recon¬ 
struction  appears  to  be  less  important,  given  that  a  voxel  size  of 
60  nm  x  60  nm  x  60  nm  is  largely  sufficient  to  describe  accurately 
the  porosity  phase.  Indeed,  in  the  studied  case,  the  microstructure 
is  relatively  coarse  in  comparison  to  the  functional  electrode.  Con¬ 
sidering  the  mean  pore  diameter  (2 Rp )  =  2.4  p,  m,  one  characteristic 
void  of  the  material  will  be  depicted  in  one  direction  by  40  elements 
of  60  nm  or  1 0  elements  of 240  nm  (leading  respectively  to  ~33,51 0 
and  ~523  voxels  included  in  the  pore  volume).  Therefore,  both  the 
RVE  size  and  the  resolution  of  the  reconstruction  are  consistent 
with  the  characteristic  length  of  the  microstructure.  As  a  conse¬ 
quence,  the  methodology  based  on  hard  X-ray  nano-tomography 
experiments  appears  to  be  well  adapted  to  investigate  SOFCs  sup¬ 
port.  Indeed,  this  method  allows  investigating  large  volumes  at 
intermediate  resolution.  Conversely,  the  3D  reconstructions  based 
on  the  FIB-SEM  technique  would  not  be  suitable,  because  (i)  the 
examination  of  large  volumes  is  time-expensive  and  (ii)  the  high 
resolution  provided  by  this  technique  is  not  necessary.  As  a  con¬ 
clusion,  it  is  proposed  to  use  the  hard  X-ray  nano-tomography 
experiments  to  obtain  reconstruction  from  SOFCs  support;  whereas 
the  FIB-SEM  technique  can  be  restricted  to  obtain  reconstruction 
of  the  functional  electrodes. 

Tortuosity  factor  -  It  has  been  proposed  in  Section  3.3  to  deter¬ 
mine  the  value  of  the  tortuosity  factor  with  the  intermediate 
boundary  conditions  depicted  in  Table  2:  this  condition  corre¬ 
sponds  to  one  flux  and  one  concentration  imposed  as  boundary 
conditions  and  leads  to  rz  =  2.8+g This  value  is  bounded  by  an 
upper  and  lower  values  corresponding  to  the  first  and  the  third 
boundary  conditions  given  in  Table  2  (i.e.  2  concentrations  or  2 
fluxes  imposed).  In  other  word,  the  tortuosity  factor  is  determined 
with  an  error  of  +11%/-  7%.  Therefore,  although  it  is  not  discussed 
elsewhere  for  SOFC  application,  it  is  worth  underlining  that  the 
calculated  effective  properties  of  porous  electrodes  has  to  be  deter¬ 
mined  by  estimating  the  bias  introduced  by  the  choice  of  the 
boundary  conditions. 

It  can  be  noticed  that  the  tortuosity  factors  determined  in  the 
present  study  are  slightly  higher  than  the  values  reported  in  litera¬ 
ture  for  Ni-YSZ  electrodes.  For  instance,  Wilson  et  al.  [17]  and  Iwai 
et  al.  [20]  have  obtained  respectively  (rx,  ry,  rz)  =  (2.1 , 2.2, 1 .9)  and 
(r*,  t y,  rz)  =  (2.03,  2.06, 1.83);  whereas  Grew  et  al.  [26]  have  found 
(r*,  t y,  rz)  =  (1.77, 1.51,  -).  These  discrepancies  in  the  results  can  be 
first  of  all  attributed  to  the  difference  in  microstructures  and  layer 
thicknesses  between  the  cermet  used  as  a  “functional”  electrodes 
or  a  cell  support.  It  can  be  also  mentioned  that,  in  our  case,  the  “non¬ 
steady”  tortuosity  factor  decreases  with  decreasing  the  size  of  the 
analysed  volume  (for  volume  lower  than  the  RVE)  before  reaching 
the  steady  value  (Fig.  8).  For  example,  an  “non-steady”  tortuosity 
factor  of  rz  =  2.01  is  obtained  on  a  volume  of  the  reconstruction 
limited  at  12  pan  x  12  pan  x  12  pm. 

Isotropy  of  the  microstructure  -  As  shown  in  Fig.  4,  the  covari- 
ograms  I<(h)  computed  in  the  three  directions  of  the  space  are  very 
similar.  That  means  that  the  characteristic  dimensions  of  the  pores 


are  equal  whatever  the  considered  direction  and  suggests  that  the 
shape  of  the  voids  is  isotropic.  This  result  justifies  the  assumption 
taken  in  Section  3.2  to  describe  the  gas  phase  by  an  assemblage  of 
spherical  pores.  Nevertheless,  the  slight  difference  in  tortuosity  fac¬ 
tors  (see  Table  2)  could  indicate  that  the  medium  exhibits  a  slight 
anisotropy  in  the  percolated  pathways  for  the  gas  diffusion. 

If  we  consider  now  the  medium  composed  by  an  assemblage  of 
mono-sized  solid  spheres,  the  tortuosity  factor  r  for  the  gas  phase 
is  given  by  the  following  expression  [42]: 


The  use  of  this  equation  leads  to  a  tortuosity  factor  of  only 
r  =  1.26  for  £  =  0.47.  It  can  be  mentioned  that  this  value  is  in  good 
agreement  with  the  ones  simulated  by  Zalc  et  al.  [43]  for  spheres 
packing  with  the  same  void  fraction.  However,  this  value  is  much 
lower  than  the  one  measured  in  this  work.  It  is  suggested  that 
the  topology  of  the  solid  phase  is  poorly  described  by  a  packing 
of  solid  spherical  particles  for  the  studied  cermet.  Rather,  solid 
phase  in  Ni-YSZ  cermets  may  be  better  represented  by  a  continu¬ 
ous  “skeleton”  embedding  spherical  pores.  In  such  pores  network, 
some  constriction  effects  may  rise  the  value  of  the  tortuosity  [43]. 

5.  Conclusion 

A  methodology  has  been  proposed  in  this  work  to  obtain  3D 
reconstruction  of  porous  SOFC  supports.  Synchrotron  X-ray  nano¬ 
tomography  measurements  have  been  carried  out  on  commercial 
Ni-8YSZ  cermet.  The  beam  energy  has  been  selected  at  29  keV  in 
order  to  obtain  simultaneously  (i)  the  required  high  transmitted 
X-ray  flux  and  (ii)  the  optimised  solid/void  phase  retrieval  condi¬ 
tions.  Thanks  to  this  experimental  procedure,  a  large  volume  of 
the  porous  material  has  been  inspected.  A  3D  reconstruction  of 
42  pan  x  42  p,m  x  105  pm  has  been  obtained  with  a  voxel  resolu¬ 
tion  of  60  nm. 

The  3D  reconstruction  has  been  numerically  analysed.  A  good 
isotropy  of  the  microstructure  has  been  highlighted  by  calculating 
the  covariograms  in  the  three  directions  of  the  reconstruction.  Con¬ 
sidering  the  morphological  properties  controlling  the  gas  diffusion, 
it  has  been  shown  that  a  volume  as  large  as  35  pm  x  35  pun  x  35  pun 
is  required  to  be  statistically  representative  of  the  whole  support. 
A  porosity  of  £  =  47%  and  a  mean  pore  radius  of  {Rp )  =  1 .2  p,  m  have 
been  found  for  the  studied  Ni-8YSZ  cermet.  The  tortuosity  factor  r 
has  been  also  calculated  on  the  basis  of  finite  element  simulations 
performed  on  the  3D  reconstruction.  It  has  been  shown  that,  even 
for  a  computed  volume  of  36  pm  x  36  pun  x  36  pun,  the  value  for 
the  tortuosity  factor  is  not  completely  independent  of  the  bound¬ 
ary  conditions  chosen  for  the  simulations.  In  this  condition,  it  is 
proposed  to  determine  the  tortuosity  factor  for  an  intermediate 
boundary  condition  which  is  relevant  for  SOFC  application.  The 
lower  and  upper  limit  of  this  parameter  can  also  be  specified  by 
considering  two  uniform  concentrations  or  flux  as  boundaries  con¬ 
ditions.  By  using  this  approach,  the  tortuosity  factor  of  the  studied 
cermet  has  been  determined  in  the  direction  perpendicular  to  the 
anode/electrolyte  interface  to  rz  =  2. 8+0  2. 
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